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ABSTRACT 

We consider the survivability of planetary systems in the globular cluster 47 Tucanae. 
We compute the cross sections for the breakup of planetary systems via encounters 
with single stars and binaries. We also compute the cross sections to leave planets on 
eccentric orbits. We find that wider planetary systems (d ^ 0.3 AU) are likely to be 
broken up in the central regions of 47 Tucanae (within the half-mass radius of the 
cluster). However tighter systems, and those in less-dense regions may survive. Tight 
systems will certainly survive in less-dense clusters where subsequent surveys should 
be conducted. 

Key words: stellar: evolution - globular clusters: general. 



1 INTRODUCTION 

The search for extra-solar planets came to fruition in 1992 
when Alex Wolszczan (1992) discovered the first confirmed 
extra-solar planets orbiting the pulsar 1257+12. Since Wol- 
szczan's discovery, conventional optical searches have found 
a number of extra-solar planets, bias ed towards the so-called 
" hot Jovians", typifie d by 51 Peg (Mayor & Queloz 1995 



Marcy & Butler 1996). The observing techniques, measure- 
ment of low amplitude periodic Doppler shifts in the pri- 
mary, bias such searches towards high mass planets with 
short orbital periods. Recently the nature of these systems, 
gaseous giant planets in short period orbits, was confirmed 
by the observation of a transit in the HD 209458 system 
( Charbonneau et al 2000 ). 

A second pulsar-planet system has also been observed. 
PSR B1620-26 in the globular cluster M4 is a binary pulsar 
with a low mass white dwar f secondary in a ~ 1 AU near cir- 
cular orbit (Lyne et al 1988). In 1993 Backer et al discovered 
the presence of a second or biting companion, conjectured 



199S 



to b e have a planet mass (Backer et al 1993; Sigurdsson 



Thorsett et al 1993). Subsequent observations have 



confirmed the presence of a second companion around 1620- 
26. The second companion has sub- stellar mass and it is near 
certain that the object is a planet (Thorsett et al 1999). M4 
is a modest, relatively low metallicity cluster, and the most 
plausible formation scenario requires that the planet was 
exchan ged into the pulsar system during an encounter with 



a main ^cqMnc^^^^TT^^l^^W^igu^Sso^^^9J ; Sig- 



urdsson 1995; Joshi fc Rasio 1997] ; 



Ford et al 



200(H ). 



The 

implication is that Jovian mass planets must exist around 
at least some of the main-sequence stars in relatively low 
metallicity globular clusters. The progenitor system of the 
M4 planet must have had the planet in a wide orbit just 



before the exchange encounter (a ^> 1AU), in order for it 
to have been plausibly exchanged into the current configu- 
ration. 

Recently, a Hubble Space Telescope project (GO-8267) 
carried out a transit search for "hot Jovians" in the globular 
cluster 47 Tucanae. The observations involved the continu- 
ous monitoring of several tens of thousands of stars using 
WF/PC2 and STIS. The basic technique was to take short 
(few minute) exposures, primarily alternating I and V band. 
A single field in the cluster, with the PC chip centered on 
the edge of the cluster core, was monitored continuously for 
over 8 days. 

The search is sensitive to Jovian sized planets with or- 
bital periods ^ 4 days, orbiting main-sequence stars in the 
cluster. A detection requires two observations of a dimming 
of a few millimagnitudes in both bands, lasting several hours 
(typical transit times scales are 3 hours). Only relative pho- 
tometry is required, and about 34,000 stars were monitored 
simultaneously. 

If the distribution of orbital parameters of Jovian plan- 
ets in 47 Tuc is comparable to that in the stellar neighbour- 
hood, for short orbital periods, then given a random orienta- 
tion of orbital planes, the search should discover some 15-20 
51 Peg like planets. 

Globular clusters are crowded environments. Binary 
systems undergo dynamical interactions with nearby stars, 
and close encounters can lead to large changes in orbital 
parameters, including the disruption of the system (Heggie 



1975 



Hut & Bahcall 198S ) . It is of interest to consider the 



probability of a planetary system surviving for a time com- 
parable to the Hubble time in the inner regions of globular 
cluster, and the possible dynamical evolution of high mass 
ratio binaries in such environments. Considerable work has 
been done on this problem, including the pioneering analy- 
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sis of Heggie (1975) and the numerical experiments of Hills 
(1975) and Hut & Bahcall (1983), and more recently sim- 
ulations directed at specific problems in globular clusters 
flDavics fc Bcnz 19951 [P 

avies 1995; pigurdsson fc Phinneyj 
1993 ); |5|igurdsson fc Phiuncy 1995^ , Relatively little work 
has been done on the problem of very high mass ratio sys- 
tems appropriate to planetary encounters, Sigurdsson (1992) 
considered the basic problem, Laughlin & Adams calculated 
mean disruption rates for open clusters (1998), and Bonnell 
& Kroupa (1998) looked at the issue in the context of planet 
formation. Hills (1984) considered close encounters between 
a star-planet system and an intruder star. Hills & Dissly 
(1989) extended this work to consider a range of impact ve- 
locities and intruder masses. For the M4 system, consider- 
able effort has been made to simulate its possible dynamical 
histo ry and the impli c ations of the exi s tence of such a sys- 
tem (sigurdsson 1993 ; Sigurdsson 1995; Joshi & Rasio 1997 
Ford ct al 2000| ). 

47 Tucanae is a massive and quite dense cluster. The 
central density is about 1.5 x 10 5 pc -3 with central veloc- 
ity dispersion of ~ 12kms _1 . The core radius is about 0.5 



pc and the half-mass radius a bout 3.9 pc (Djorgovski 1993 
Howell & Guhathakurta 2000). At the half mass radius the 
number density of stars is about a factor of ten lower than 
in the core. 



2 ENCOUNTER TIMESCALES 

Encounters between two stars will be extremely rare in the 
low-density environment of the solar neighbourhood. How- 
ever, in the cores of globular clusters number densities are 
sufficiently high (~ 10 5 stars/pc 3 in some systems) that en- 
counter timescales can be comparable, or even less than, the 
age of the universe. In other words, a large fraction of the 
stars in these systems will have suffered from at least one 
close encounter or collision in their lifetime. 

The cross section for two stars, having a relative velocity 
at infinity of Voo , to pass within a distance R m in is given by 



VI 



(i) 



where V is the relative velocity of the two stars at closest 
approach in a parabolic encounter (i.e. V 2 = 2G(M\ + 
M 2 )/R m i n ). The second term is due to the attractive grav- 
itational force, and is referred to as gravitational focussing. 
If V 2> Voo as will be the case in systems with low velocity 
dispersions, such as globular clusters, a oc -R m i n . 

One may estimate the timescale for a given star to un- 
dergo an encounter, r cnc = l/nav. For clusters with low 
velocity dispersions, we thus obtain 



7xl0 10 yr 
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where n is the number density of single stars of mass M. 

We may estimate the timescale for an encounter be- 
tween a star-planet system and a single star, in a similar 
mariner where now i? m i n — d, where d is the semi-major axis 
of the planet's orbit. The encounter timescale for such a sys- 
tem may therefore be relatively short as the semi-major axis 
can greatly exceed stellar radii. For example, a planetary 
system with d ~ 1AU(= 216i?o), will have an encounter 



timescale r cnc <J 10 years in the core of a dense globular 
cluster. 

Encounters between single stars and binary systems will 
lead to the break up of the binaries when the kinetic energy 
of the incoming star exceeds the binding energy of the bi- 
nary. This will be the case when Voo > V c , where 

V 2 GM 1 M 2 {M 1 + M 2 + Ms) ,„v 
M 3 (M 1 +M 2 )d ( ' 

where Mi is the mass of the primary, M 2 the mass of the 
secondary and M3 is the mass of the incoming, single star. 
Binaries vulnerable to break up are often re ferred to as be- 
ing soft ( Heggie 1975 ; Hut fc Bahcall 1983 ). In the case of 
encounters between single stars and a star-planet system, 
the latter will be invariably soft as M 2 <C Mi and M3. We 
would therefore expect the planetary system to be broken up 
with some high probability, when an intruding star passed 
within a distance d of the planetary-system star (but see 
Sigurdsson 1993). In addition, more distant encounters may 
leave the planets in somewhat eccentric orbits, as will be 
discussed later. The importance of encounters between a 
planetary system and binaries is a function of the binary 
fraction within the cluster, which is not well known. 



3 NUMERICAL METHODS 

For encounters involving a single intruding star, we used a 
code based on the one discussed in Davies, Benz & Hills 
(1993). For encounters involving binaries (i.e. four objects), 
we used a code based on the one discussed in Bacon, Sig- 
urdsson & Davies (1996). 

The initial conditions for the scatterings were set follow- 
ing the method of Hut & Bahcall (1983, see also Sigurdsson 
& Phinney 1993). With two binaries, we have additional 
parameters from the relative phase of the second binary, 
the orientation of the plane of the second binary and the 
second binary mass ratio, semi-major axis and eccentricity. 
We drew the binary parameters by Monte Carlo selection 
uniformly over the phase variables. The relative velocity at 
infinity of the centres-of-mass of the two binaries, Uoo was 
chosen uniformly on the interval allowed. The initial binary 
eccentricities, ei,2 were zero for all encounters, previous cal- 
culations indicate the cross-sections of interest are not sen- 
sitive to the binaries' eccentricities. 

For encounters involving single stars, the equations of 
motion are rescaled using the Heggie (1972) quasi regulari- 
sation and are integrated using the variable-order, variable- 
stepsize Shampine-Gordon integrator (1975). For encounters 
involving binaries, a Bulirsch-Stoer variable step integrator 
with KS-chain regularisation (Aarseth 1984, Mikkola 1983, 
1984a,b) was used for integrating the motion of the particles. 



4 RESULTS 

We considered encounters between a star-planet system and 
either a single star or a binary. The star-planet system con- 
tained a single planet in a circular orbit. We consider en- 
counters for various values of a, where a = Voo/V or b. These 
will correspond to different values of Voo as a function of the 
radius of the planetary orbit, d, and thus V or b. The encoun- 
ters considered are listed in Tables 1 and 2. In both tables, 
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Mo /Mn 


CK 


v orb 


Vac 


V c 


0.1 


1.0 


31.24 


31.24 


13.01 


0.01 


1.0 


29.94 


29.94 


4.20 


0.001 


1.0 


29.80 


29.80 


1.33 


0.01 


3.0 


29.94 


89.82 


4.20 


0.01 


0.55 


29.94 


16.46 


4.20 


0.01 


0.30 


29.94 


8.98 


4.20 


0.01 


0.17 


29.94 


5.09 


4.20 


0.01 


0.10 


29.94 


2.99 


4.20 


0.001 


0.30 


29.80 


8.94 


1.33 



Table 1. A list of the binary-single star encounters studied, a = 
Voo/Vorb and V c is the critical velocity given by equation (3). M 2 
is the planet mass. All velocities are given in km/s. 



M 2 /M a V orb Voo V c 



0.01 


3.0 


29.94 


89.82 


36.54 


0.01 


1.0 


29.94 


29.94 


36.54 


0.01 


0.30 


29.94 


8.98 


36.54 


0.01 


0.10 


29.94 


2.99 


36.54 



Table 2. A list of the binary-binary encounters studied, a = 
Voo/Vorb and V c is the critical velocity given by equation (4). M 2 
is the planet mass. All velocities are given in km/s. 

the values of V OT b and Voo are given in km/s assuming d = 1 
AU. In Table 1, V c is given by equation (3) whereas in Table 
2, the velocity at which the total energy is zero, V c is given 
by 

= G J + MiM±) (4) 

\ di d 2 J 

where n = M T /(Mi + M 2 )(M 3 + M 4 ), M T = Mi + M 2 + 
M 3 + M 4 , is the binaries reduced mass, and di,d 2 are the 
semi-major axis of the binaries containing masses Mj.,2, -^3,4 
respectively. 

Note for all simulations, Mi = M 3 = M 4 = 1.0M©, and 
di — di. The turnoff mass for 47 Tuc is slightly lower than 
1 M©, probably 0.85±0.05M© , and the true encounter rates 
scale to be correspondingly slightly lower. 

For each value of a, we simulated encounters for a range 
of impact parameters, ensuring that we considered sufficient 
large impact parameters to include all relevant interactions. 
For each set of encounters we then computed the cross sec- 
tions for the planetary system to be disrupted (i.e. the planet 
does not remain in orbit around either star), and for the 
planet's orbit to be perturbed to eccentricities, e > 0.03, 
0.1, and 0.3. 

The cross section for the breakup of the planetary sys- 
tem via encounters with a single star is shown in Figure 1 
as a function of a. Also shown is the cross section for the 
two stars to pass within a distance, i? m in < d which was 
calculated using equation (1). Also shown is the breakup 
cross section for a — 0.3 rescaled assuming that the change 
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1 1 ' ' ' ' 

0.1 1 10 

a 

Figure 1. The cross section for breakup of the planet-star system 
as a function of a = Voo/V OIO (solid line), and for the intruding 
star to pass within the orbital separation of the planet-star system 
(dotted line). Both cross sections arc given in units of the binary 
separation (ie the cross section of a target having a radius equal to 
the binary separation is n). The dashed line is the breakup cross 
section for a = 0.3 rescaled assuming that the change in cross 
section is purely an effect of the differing degree of gravitational 
focussing ocurring for different values of a. 



in cross section is purely an effect of the differing degree of 
gravitational focussing ocurring for different values of a. 

The cross section for the breakup of the planetary sys- 
tem via encounters with binary stars shows a similar form 
to that seen in Figure 1, except the value is increased by a 
factor ~ 2. 

We consider now the cross sections for the production 
of eccentric planetary systems. These are shown in Figure 
2 for encounters with single stars. The cross sections are 
larger than the equivalent breakup cross section. The cross 
section for producing a planetary system with eccentricity 
e > 0.03 being ~ 3-5 Obrcakup ( see a l so Heggie & Rasio 
1996). An important point to note here is that the systems 
with non-zero eccentricities can have semi-major axes larger 
than the initial circular orbit of the planet, and in any case 
the maximum separation (i.e. d nC w(l + e)) is certainly larger. 
A planet left in such an eccentric orbit is therefore a larger 
target for a subsequent breakup. We estimate that this will 
double the effective breakup cross section for planets with 
initial separations d ~ 0.1 AU. Systems left in extremely 
eccentric orbits tend to have equivalently large semi-major 
axes, this being the limit of sytems that just failed to be 
broken up. The number of eccentric systems where circular- 
isation occurs will be small and limited to systems where 
the initial separation d <5 0.1 AU. 

If we consider that V4o is a constant value for all sets 
of encounters, then the various values of a correspond to 
different values of orbital separation d. As a = Voo/Vorb 
and Uorb oc \fljd, d oc a 2 . Considering Uoo = 9 km/s, 
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Figure 2. The cross section for the planet to be left in an eccen- 
tric orbit as a function of a = Voo/V or b. Cross sections are for 
systems with eccentricity e > 0.03, 0.1 and 0.3 (top to bottom). 
All cross sections are given in units of the binary separation as in 
Figure 1. The dashed line is the breakup cross section for a = 0.3 
rescaled assuming that the change in cross section is purely an 
effect of the differing degree of gravitational focussing ocurring 
for different values of ex. 

inspection of Table 1 reveals that a — 0.3 corresponds to 
d — 1 AU. rescaling the other cross sections for the same 
value of Voo we produce Figures 3 and 4, where now the 
cross sections are given in units of (AU) 2 . These may be 
converted to timescales via 

. = 4.16xlO-yr(I)(i) (5) 

where n is the number density of stars (or binaries) /pc 3 , and 
a is the cross sections in units of (AU) 2 . In Table 3 we give 
the timescales in years for systems to be broken up or left 
in eccentric orbits for both binary-single and binary-binary 
encounters, assuming a number density, n = 10 5 stars/pc 3 , 
and Voa = 9 km/s. We will consider the implications for the 
survival, or otherwise, of planetary systems in stellar clusters 
in the discussion section. 



Figure 3. The cross section for breakup of the planet-star system 
from encounters with single stars as a function of separation in 
AU (solid line), and for the planet to be left in an eccentric orbit 
with eccentricity, e > 0.03, 0.1, 0.3 (dashed lines, from top to 
bottom). 



10° F 



10 5 r 
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4.1 The importance of collisions 

Thus far we have only considered exchanges and breakups 
as mechanisms to destroy planetary systems. For the very 
tightest systems collisions between the planet and a star or 
between two stars will also play a role. If two stars collide 
and merge, the planet will likely be swallowed up by the 
expanding merged object, assuming it remainds bound to 
the collision product in the first place. 

We computed the cross sections for physical collisions or 
very close encounters between objects during an interaction. 
For encounters involving a single star, we find that mergers 
are unimportant except for encounters involving very tight 
planetary systems (i.e. d ^ 0.1 AU) where a co \ ~ 0.5abrcaku P . 



l 1 1 — 1 — 1 — 

0.1 1 10 100 

binary separation (AU) 

Figure 4. The cross section for breakup of the planet-star system 
from encounters with binary stars as a function of separation in 
AU (solid line), and for the planet to be left in an eccentric orbit 
with eccentricity, e > 0.03, 0.1, 0.3 (dashed lines, from top to 
bottom). 
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d/AU r breakup e > 0.3 e > 0.1 e > 0.03 



binary— single encounters 



1.0000 
0.1000 
0.0500 
0.0300 



5.2E+08 2.0E- 

5.2E+09 2.0E- 

1.0E+10 4.0E- 

1.7E+10 6.7E- 



•08 
09 
09 
09 



binary— binary encounters 



1.0000 3.0E+08 1.1E+08 

0.1000 3.0E+09 1.1E+09 

0.0500 6.0E+09 2.7E+09 

0.0300 1.0E+10 3.8E+09 



1.2E+08 
1.2E+09 
2.4E+09 
4.0E+09 



5.6E+07 
5.6E+08 
1.1E+09 
1.9E+09 



8.7E+07 
8.7E+08 
1.7E+09 
2.9E+09 



3.7E+07 
3.7E+08 
7.4E+08 
1.2E+09 



Table 3. Timescales (in years) for systems to be broken up or 
left in eccentric orbits for both binary-single and binary-binary 
encounters, assuming a number density, n = 10 5 stars/pc 3 , and 
Voo = 9 km/s. 



For encounters involving binary stars, we find that the cross 
section for mergers can be larger than the breakup cross 
section for d ^ 0.5 AU. For such tight planetary systems, 
we find that <j co i ~ 150 in units of (AU) 2 . 



5 DISCUSSION 

Combining the cross sections for breakups, significantly in- 
creasing system size (by making it eccentric) , and collisions, 
we estimate that the effective cross section for the destruc- 



30 for encounters 
150 for encounters in- 



tion of tight planetary systems (Xdestr 
involving single stars, and (Xdestr 
volving binary stars. Assuming the core binary fraction in 
47 Tuc to be 0.1 - 0.2, we therefore conclude that effective 
destruction cross section in the core is ^destroy ~ 40 — 50, for 
a planet in 0.4 AU orbit. The data from the HST project will 
give the fraction of high mass-ratio main sequence binaries 
in the same part of the cluster surveyed for planets. 

Using equation (5), we are able to compute the likely 
timescale for the destruction of a planetary system, or the 
time required to place the planet on an eccentric orbit. 
Assuming that a particular star-planet system is located 
within the core of 47 Tuc, where n ~ 10 5 stars/pc 3 , we see 
that Tbreakup ~ 4 x 10 10 /^destroy years. We therefore require 
^destroy <D 4 which is easily satisfied for most semi-major 
axes. A planet with an orbital separation d ~ 5 AU would 
be broken up in ~ 10 s years. By comparison, a planet orbit- 
ing a star at the cluster half-mass radius will be somewhat 
less vulnerable. Here the number density n ~ 10 4 stars/pc 3 
hence we require <7d cs troy Jj 40; only those systems where 
d <; 0.3 AU are likely to be destroyed. The critical value 
of d will clearly depend on the exact value of the number 
density of stars, n. Note that Sigurdsson (1992) assumed a 
stellar density ~ 3 times smaller for 47 Tuc and derived a 
correspondingly longer time scale for breakup in the core. 
Additional uncertainty is introduced by the time evolution 
of the 47 Tuc stellar density profile. If the cluster is evolv- 
ing through normal relaxation, then it was substantially less 
dense in the centre a few billion years ago. However, there is 
some suggestion that the cluster may already have evolved 



through a dense core phase ( {Sills et al 2000| ) in which case 
rate of destruction of tightly bound planetary systems may 
have been an order of magnitude higher in the recent past. In 
either case, theory predicts little or no evolution in the stel- 
lar density near the half-mass radius, and the WF/PC2 field 
extends out from the core out almost to the half-mass ra- 
dius. On the other hand, there are more stars to be oberved 
in the denser parts of the cluster and the prior expectation 
therefore that most planets should have been found in the 
densest part of the observed field, assuming no dynamical 
disruption took place. 

It should also be noted that simulations of the evolution 
of multiple planetary systems suggest that a relatively slight 
perturbation in the eccentricity of one or more planets may 
lead to a radical rearrangement of the system i n a short 
timescale ( Quintan 1992 ; Murray fc Holman 1999 ) Thus we 
note that encounters that perturb the planets may play a 
role in the evolution of planetary systems, even in relatively 
low-density short-lived stellar clusters. For example, if n ~ 
10 3 stars/pc 3 , r cnc ~ 10 7 years, for d ^ 30 AU. 

In terms of the current survey in 47 Tuc, it would seem 
possible that all planetary systems have been destroyed if all 
the stars in the survey were located in the highest-density 
regions. Given that some of the stars observed are in regions 
with an order of magnitude smaller stellar density, it is dif- 
ficult to explain away all the expected detections. However, 
there are fewer stars observed in the low-density regions of 
the cluster, so if the intrinsic specific density of 51 Peg like 
systems is somewhat lower than in the solar neighbourhood, 
then we might have expected more like 5-10 detections from 
observing 34,000 stars, and 75% of those in the densest re- 
gions. In order for direct dynamical disruption to explain 
the absence of 51 Peg like systems, one still requires either a 
past denser phase for the core, or a very long time spent at 
the current density. The minimum average density for the 
environment of a given planetary system required for its de- 
struction being ~ 10 4 stars/pc 3 . From our calculations here, 
we conclude that tight planetary systems (i.e. d & 0.1 AU) 
will definitely survive in less-dense clusters where subsequent 
surveys should be conducted. An ideal target is NGC 6352. 
It has core density of about 10 3 stars/pc 3 , is at an estimated 
distance of 6.1 kpc (compared to 4.6 kpc for 47 Tuc) and has 
metallicity [Fe/H] = -0.51 (compared with -0.71 for 47 Tuc). 
With the Advanced Camera on the Hubble Space Telescope, 
a considerably wider field can be monitored, and another 
survey for transits carried out. 
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